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1 

The effect of disorder on a class of transition metal oxides described by a single orbital Hubbard 

^ ' model at half filling is investigated. The phases are characterized by the nature of the electronic and 

(N 
C"^^ . spin excitations. The frequency and temperature-dependent conductivity and spin susceptibility 

m 

^^ I as functions of disorder are calculated. The interplay of disorder and electron-electron interaction 

(^ ' produces unusual behavior in this system. For example, the dc conductivity, which is vanishingly 

OO ■ 

^D ■ small at low disorder in the Mott phase and at high disorder in the localized phase, gets surprisingly 



enhanced at intermediate disorder in a "metallic" phase. Moreover, the spin susceptibility in this 



H ' "metallic" phase is not the expected Pauli-behavior but Curie-l/T due to the presence of local 

moments. 

PACS numbers: 



I. INTRODUCTION 

Based on the scaling theory, it is expected that non-interacting electrons are localized for 
any weak disorder in two dimensions in zero magnetic field [l|. This theoretical result was 
questioned when transport experiments in 2D silicon metal-oxide-semiconductor field-effect 
transistors (MOSFET's) showed a metallic phase with ^ > 0, where p is the resistivity, at 
higher density of carriers than the critical density of carriers, ric and an insulating phase 



with -1^ < at lower density of carriers, n < Uc 



. Further experiments supported the 



observation of a MIT in 2D systems [d, m, U U U Isl, bl, M- 

The interplay of disorder and electron-electron interaction can have an important effect 
in 2D systems. The recent experiments have shown that electron-electron interactions can 
lead to a possible metallic phase in the 2D electron gas. Motivated by these experiments, 
we investigate the 2D disordered Hubbard model at half- filling pjj]. In previous work, a 
possible metallic phase at T = with extended wave functions, sandwiched between a 
Mott insulator for low disorder and a localized insulator for large disorder, is found. In the 
metallic phase, the spectral gap is zero but antiferromagnetism is found to persist. The 
existence of a metallic phase is confirmed by the finite size scaling behavior of the inverse 
participation ratio (IPR). At intermediate disorder, the IPR for f/ = extrapolates to 
a finite value as L ^ oo. Thus, the system is a localized. However, for the interacting 
system, U = At, the IPR as a function of 1/L extrapolates to zero in the thermodynamic 
limit, implying a diverging localization length would occur in a metallic phase. The fact 
that the non-interacting system with the same disorder realization had a finite localization 
length indicates that interactions clearly played an important role in delocalizing the wave 
functions or at least in making them much more extended than the wave functions in the 
absence of interaction. 

The study of the Hubbard model has been stimulated after the discovery of cuprate su- 
perconductors, such as YBa2Cu307_x where x is the doping concentration and x = 1 — n 
where n = -^f^^ is the concentration of electrons [l2]. A systematic phase diagram of the 
copper oxide materials is shown in Fig. [1] where x = corresponds to half-filling with 1 elec- 
tron per site. At x = the system is a Mott insulator with anti-ferromagnetic order. Upon 
doping, the Neel temperature, T/v, decreases. On the other hand, T^ shows a dome shaped 
dependence on x. This doping not only gives superconductivity but also introduces disorder. 
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FIG. 1: Systematic phase diagram of the copper oxide materials. 



A realistic modeling will therefore need to consider both effects of electron correlation and 
disorder. 

Our focus is to understand the competition between random disorder and electron- 
electron interaction on the Mott phase with long-range AF order in 2D. In order to charac- 
terize the phases, we calculate the temperature and frequency dependent conductivity and 
spin susceptibility. The calculation are done for the two-dimensional disordered Hubbard 
model on lattices of size up to 32 x 32 at T = and as function of temperature and disorder 
using a numerical Hartree-Fock self-consistent method. Our results are reported at half 
filling for interaction strength, U = At. 

The remainder of this paper is organized as follows. In Section [Til we give the model and 
the Hartree-Fock method. In Section llllt the self-consistent numerical method is presented. 
In Section [IVl we give our results of the conductivity and spin susceptibility. This is followed 
by a concluding discussion and an outline of possible future research. 



II. MODEL 

We consider the disordered strongly correlated electron system described by the 2D one- 
band Hubbard model with site disorder: 

H = -tY^ [clcj^ + cJ^Q^j + Y^ Uui^Uii + ^(Ki - fi)clci^ (1) 

where t is the near-neighbor hopping amplitude, U is the on-site repulsion between electrons, 
cjcr {cia-) is the creation (destruction) operator for an electron with spin a on site r^, Ujo- = 
CifjCia is the number operator at site r^, and fi is the chemical potential. Vi is local disorder 
potential at site Tj and is chosen from a uniform distribution, [—V", V], where V is the measure 
of the strength of disorder. 

The Hubbard model with U ^ lOt describes the essential physics of the copper oxide, 
Cu02, perovskites upon doping. At half filling, (n) = 1, with no disorder, the system at 
T = has long-range antiferromagnetic (AF) order and is a Mott-type insulator. 

The mean-field approximation is applied for the electron-electron interaction with the 
following prescription: 

n^m^ ~ {ni^)n^ + {n^)nii - {ni^){ni^) - c\^Ci^{c\ci^) - c\^Ci^{c\^Ci^) + {c\^Cii) {c\ci^) . (2) 

Then the effective Hamiltonian, ignoring constant terms, — X]jf^{(^it)(^«i) + ('^lT'^«i)('^li'^«T)}' 
becomes 

{ij),cr i,<y i,cr i 

where h^ = —U{c\^Cii), h^ = — f/(c|^Ci|), f =i, and | =]. In our calculation, we assume hi 
and h~ are real so hi = hf = h~ . We therefore have 3N variational parameters, (rii^) for 
both spin species and hi at the A^ site that must be determined self-consistently. 

Exphcitly, the effective Hamiltonian for A^ lattice sites in the bases, llf), |1|), |2t), |2|) ... 



lA^f)) l^i) with periodic boundary conditions is: 
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Vi + U{ni^) - fi. 
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III. METHOD 



The inhomogeneous Hartree-Fock approximation for the disordered Hubbard model gen- 
erates a matrix with local densities, {ni^j), and local magnetic fields, hi, which are variational 
parameters and must be solved self-consistently for each /i |llUl3l|. For half filling, /i = U/2 
with no disorder but in the presence of disorder, fi is adjusted to satisfy the half filling 
condition. The initial input parameters for {riicr) and hi are determined as follows: For low 
disorder we start the self-consistent procedure with an AF initial condition for (riicr). On 
the other hand, for high disorder a paramagnetic state is used for the stating configuration. 
For both cases, hi are chosen from an AF initial condition from — O.lt or O.lt. The initial 
condition for hi should not affect the results. The Broyden method is used for achieving self- 
consistency efficiently [ij]. Input and output fields are compared after each iteration and if 
the difference of the fields at all sites is less than 10~^, the self-consistentcy loop is exited. 
Using these final self-consistent fields, {riia) and hi, and eigenvalues, physical quantities of 
interest such as conductivity are calculated. We also test that the self-consistent values of 
the parameters are not dependent on the initial starting conditions for the parameters. 



IV. RESULTS 

Our main effort in this part of the thesis is to calculate the frequency and temperature 
dependent conductivity and spin susceptibility, which are experimentally measurable and 
provide further characterization of these phases. 

Summary of results are the following. In Region 2, the conductivity near cu = is most 
enhanced and shows the Drude like behavior and the peak frequency, where the conductivity 
has the maximum value, is lowest. The Drude model can explain non-monotonic behavior of 
the dc conductivity as a function of T. The result of ac conductivity shows non- Fermi liquid 
type behavior in Region 2. The sum rule is justified except for high T at low disorder. The 
temperature dependent behavior of spin susceptibility, Xi shows the nature of moments. At 
high T, the moments are free as indicated by their Curie behavior; at low T, x is strongly 
suppressed because electrons get paired by forming singlet on low disorder and by on-site 
pairing on high disorder sites. 

A. Frequency dependent conductivity at T = 

The real part of frequency dependent conductivity, o"(ct;), is related to the imaginary part 

I I 
of the current-current correlation function [l^, 

Re|.M| = !=^^. (5) 

where A(c(j) is the Fourier transform of A(r) = {j{t)j{0)) and j is the paramagnetic current 
density operator. Fig. [2] shows cr{uj) for various disorder strengths and Fig. [3] shows the gap 
of the conductivity as a function of disorder and the value of u where the conductivity has a 
peak as a function of disorder. For low disorder. Region 1, there is a gap in the conductivity 
which is twice the spectral gap seen in the DOS and a strong peak in the absorption at 
u c^ U. At intermediate disorder, the gap in a^uj) closes at u; = in Fig. [3] (a) but the peak 
in Re[cr(co')] persists. As seen in Fig. |3] (b), the peak moves to lower values of u as disorder 
increases in Region 2. In Region 3, the peak starts moving out to higher frequencies but 
the conductivity at high T is still steeper than Re[cr(co')] ~ l/tu^. With increasing disorder, 
the spectral weight above the Mott gap is transfered to lower spectral region, enhancing 
the absorption at a low frequency. At Vd = 2.2t, the gap of the conductivity closes because 
energy scale of disorder is comparable with the gap. Near Vd = 3t, Re[cr(u; = 0)] is most 
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FIG. 2: The frequency dependence of the conductivity, cr{uj), for V^ = 1, 2, 3, 4, and 5. The system 
size is 32 X 32, [/ = 4f , and the data is averaged over 20 disorder reahzations. The low frequency 
region of cr(u;) for V^ = 3t is obtained by a third order polynomial fit of the current-current 
correlation function, A = uja. For Vd = It there is clearly a finite Mott gap at the low frequency. 
Thus, cr{uj) is suppressed at low frequency and there is a peak in the frequency when to is around 
the Mott gap. With increasing disorder the spectral weight above the Mott gap is transfered to 
lower spectral region and the peak of cr(uj) moves to the low frequency. At V^ ~ 2.2t, the gap in 
the conductivity closes because energy scale of disorder is comparable with the gap. Near Vd = 3t, 
the conductivity near u; = is most enhanced and the frequency where the conductivity has the 
maximum value is lowest. For Vd = 4t and Vd = 5t, the dc conductivity, cr(0), is zero due to the 
Anderson insulating regime. In Region 2, the peak location of cr(u;) moves to the low frequency 
as the Mott gap closes but with increasing disorder the peak location again moves to the high 
frequency in Region 3 



enhanced and ujpeak, frequency where the conductivity has the maximum value, is lowest. 
The optical conductivity at V^ = 3t shows a Drude-like behavior. As disorder increases, this 
Drude-like behavior disappears. The origin of this Drude-like behavior at V^ = 3t may be 
the result of the screening of the strongly disordered paramagnetic sites with 2 electrons on 
a site. The disordered potentials are screened and an electron feels the following density 
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FIG. 3: (a) The gap in Re[(T(aj)] as a function of disorder, Vd- (b) ujpeak as a function of disorder. 
The system size is 32 x 32, U = 4t, and the data is averaged over 20 disorder reahzations. tOpeak is 
obtained by fitting Re[cr(u;)] in Fig. linear the peak using Lorentzian. The gap in the conductivity 
clearly closes at V^ ~ 2.2t. ojpcak in Region 2 is lower than that in Region 1 and 3. 

dependent effective potentials. 

V,^ = Vi-fi + U{n,^) 
Vil = Vi-i^ + U{ni^) 

where the effect to an electron comes from the density with opposite spin. Ignoring spin- 
polarization, this effective potential becomes 



V = li% + V,^)=V,-fi+^ ((n,t) + (n,;)) . 



(6) 



For Vd = 4:t and Vd = 5t, although the Mott gap is completely filled, the conductivity is 
suppressed due to the Anderson insulating regime. This insulating behavior continues and is 
enhanced more as disorder increases. It is difficult to fit the conductivity at high frequency 
by the form (l/u;°) and o"(u;) for all case has a much faster decay with a = 2, which is a 
conventional Fermi liquid result. 

B. DC conductivity 



DC conductivity is obtained by 



adc = lim, 



■uj^O' 



ImA(Li;) 



LU 



(7) 




FIG. 4: Temperature dependence of the dc conductivity, add^), for Vd = 1,2,3, and 4. The 
system size is 16 x 16, C/ = 4t, and the data is averaged over 20 disorder reahzations. The dc 
conductivity at low T for V^ = 3t is obtained by a third order polynomial fit of the current-current 
correlation function, A = ua. At low T, the dc conductivity in Region 1 is zero because of the 
Mott gap, the system with Vd = 3t is more conducting than Vd = It due the close of the Mott gap, 
and with increasing disorder, adc is suppressed because of an Anderson localization of electrons. 
On the other hand, at high T, adc in Region 1 is higher than the others because energy scale of 
temperature is comparable with the Mott gap and as disorder increases, adc decreases due to the 
increase of transport scattering of electrons and disorder. It is demonstrated that disorder can 
drive a metallic phase in the system. At low temperature, the dc conductivity in Region 1 and 2 
decreases as lowering temperature, signifying an insulating phase. However, the dc conductivity 
for Vd = 3t has been most enhanced and slope of the conductivity, adc/dT, has no slope. 



Fig. m shows the temperature dependence of adc for various V. In the Mott insulator (e.g. 
Vd = It or Yd = 2t) , adc is exponentially small at low T, it peaks at a temperature of order 
U , and at higher T, adc decreases with T. In the intermediate disorder regime the system 
with Vrf = 3t is more conducting than Vd = It at low T but at high T it is the opposite 
way. Increasing the site disorder further drives the system from a metallic to an Anderson 
insulating phase and the low T, adc drops again. The dc conductivity at low T and at low 



disorder e.g. < V^ < 2 is zero because of the Mott gap, adc at moderate disorder e.g. 
2 < V < 3.5, has a finite value, and then again adc is zero due to the locahzation of electron 
on the strongly disordered site. The dc conductivity shows non-nionotonic behavior as a 
function of T which can be understood using the Drude model: 

<^r'' = ^, (8) 

m 
where r is the relaxation time and n is the density of electrons. At high temperature 
adc is dominated by decreasing r. For increasing disorder, r decreases because of enhanced 
scattering between electrons and the disorder potential. At low temperature adc is dominated 
by the temperature dependence of n. The number of electrons need to be excited beyond 
the Mott gap so the density is proportional to q-^mou/t ^ Therefore, the conductivity at low 
disorder increases rapidly with increasing temperature. For both low and high disorder, adc 
vanishes as T ^ 0, signifying an insulating phase, either a Mott insulator or an Anderson 
insulator. However, as T — > the dc conductivity for V^ = 3t is unusual; it is finite with 
negligible T-dependence at low T. 

C. Frequency dependence of the conductivity for finite temperature 

Fig. O shows the frequency dependence of the conductivity, Re[cr(u;)], for various disorder 
strengths. Region 1: Mott Insulator at low disorder (V=l) shows a finit gap in the spectrum 
and a peak in Re(ci;) at a characteristic frequency uj-^^ak ~ U . With increasing T, uo^cak moves 
to lower uj and the weight at Recr(u; ^ 0) up. At T ^ t the weight at intermediate frequencies 
i^peak vanishes and gets concentrated in a Drude-like peak around u; ~ 0. Region 2: Metalic 
regime at intermediate disorder (V=2,3). 'Rea{(jj) shows that the Mott gap vanished at low 
T. As T increases weight buids up at u; ~ and has a finite intarcept, adc- At higher T, 
a{uj) &i uo K, {] develops and shows Drude-like behavior. Region 3: Insulating regime at high 
disorder (V=4), adc continues to be zero even at higher temperature. Wpeak shifts to lower 
uo with increasing T. 

For low disorder at low T there is a Mott gap but with increasing T this gap closes and 
a{uj = 0) has large weight. In Region 2, for all T, the Drude like conductivity is seen and in 
Region 3, this behavior is less obvious. AC conductivity from the Drude theory is 
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FIG. 5: The frequency dependence of the conductivity for V^ = 1 (left top), V^ = 2 (right top), 
V^ = 3 (left bottom) and V^ = 4 (right bottom) at different temperatures. The system size is 
16 X 16, [/ = 4t, and the data is averaged over 20 disorder realizations. The low frequency region 
of cr{uj) for V^ = 3t is obtained by a third order polynomial fit of the current-current correlation 
function, A = toa. For low disorder at low T there is the Mott gap for low frequency and with 
increasing T, the peak in the conductivity moves to the weight at a; = 0. With increasing disorder 
from Region 1 to Region 2, the spectral weight above the Mott gap is transferred to lower spectral 
region at low T and cr(w) shows the Drude like behavior for all T. With increasing disorder from 
Region 2 to Region 3, this Drude like behavior is less clear. The Drude like behavior shows non- 
monotonic behavior as a function of T. The conductivity weight in low frequency increases with 
increasing T up to T = 0.8 but with further increasing T it decreases. 
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where cxo = ^^|^. The real part is 

a^-'^^iuj) = -—^. (10) 



At high frequency, a^'^^'^^{uj) = ne'^/{muj'^T). It is difficult to fit the conductivity at high 
frequency by the form (l/cj") because cr{uj) for all cases has a much faster decay with a = 2. 
The sudden increase of the conductivity at low u with increasing T for V^ = 1 and V^ = 2 
may be explained by the exponential increase of the DOS at the Mott insulator as a function 
of temperature. 

D. Sum rule for the conductivity 



As a test of our numerical results, we check the consistency of the f-sum rule 16|] . The 
frequency integral of Re[o"(ci;)] is related to the total kinetic energy: 

^"^^ -A rReKu;)]rfcu, (11) 



N vre- JO 

where T = t J2{ij) a '^L'^jo-- We calculate and compare the left-hand side and the right-hand 
side of Eq. ( ITTI) . independently in Fig. [6l These values are non-monotonic as functions 
of both T and disorder. The average kinetic energy, (T), increases with temperature at 
low T but decreases with T for higher temperature. Similarly, with disorder, the kinetic 
energy increases for low disorder around the Mott insulator but at high disorder the kinetic 
energy decreases with increasing disorder because the mobility of electrons is decreased at 
strongly disordered sites. At low temperature, the sum rule is in good agreement. At finite 
temperature the agreement in the sum rule improves at higher disorder. The disagreement 
at high temperature for V^ = 1 comes from the difficulty of calculating exactly the Drude 
weight, which has a singular behavior at a; = 0. 

E. Uniform static spin susceptibility 

The frequency-dependent spin susceptibility is defined by 

r'iq,uj)= [ e^-V^(g,r)dr, (12) 

Jo 
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FIG. 6: Comparison of -^^^ (red square) where f = t ^/j,\ o- cj^Cja and ^ f^ a{uj)doj (blue dot) 
as (a) a function of disorder for T = or a function of a temperature for (b) V^ = 1, (c) V^ = 2, 
and (d) V^ = 3. The system size is 16 x 16, U = 4t, and the data is averaged over 20 disorder 
reahzations. The agreement between the blur and red dots indicates that the f-sum rule is obeyed. 
The disagreement in (b) at low disorder (V=l) comes from the difficulty of calculating the Drude 
weight in the presence of singular (5- function at u; = 0. 



where 



X^^(g,r) = (5,(g,r)5,(-g,0)) = 1 ^e^«^-(-/)(5;(r)^|(0)) 



(13) 



«j 



and Sz is z component of the spin operator. Fig. [7] shows the uniform spin susceptibihty, 
x(0, 0), and staggered spin susceptibihty, ^(vr, vr), for various disorder strengths as a function 
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FIG. 7: x(Oi 0) and xk^^ ^) as a function of a temperature for V^ = 0, 1, 2, 3, 4, and 5. The system 
size is 16 x 16, [/ = 4t, and the data is averaged over 20 disorder reahzations. The tails of x(0, 0) 
are fitted by the function ^^r^ where C is the measure of local moments and a is the measure of 
the strength of paramagnetism. These C and a are plotted as a function of the disorder, V. At low 
temperature, x(0, 0) and %('''"; ^) foi" low disorder are suppressed and these are most enhanced in 
Region 2, and with increasing disorder further, x(0,0) and x('?'";7r) decrease. On the other hand, 
at high temperature, increasing disorder, x(0, 0) and x(7i";7r) monotonically decrease. 

of T. At low temperature, for low disorder, x(0, 0) vanishes indicating a spin gap. x(7r, tt) is 
finite due to long-range AF order in the system. As T increases x(0, 0) and xk^ ■, ^) increase, 
reach a maximum at a characteristic temperature and decrease at higher temperatures. The 
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spin susceptibilities are most enhanced at intermediate disorder, (Vd = 3t). 

Also, the fitted parameters, C and a, where ^^ is used for the fitting the tail of x(0, 0), 
are considered to see local moments and the strength of paramagnetism. On the other hand, 
the Curie susceptibility of the spin at ksT ^ q^irH for localized electrons is 

Curie _ n{giiBfS{S+l) 
^ 3kBT ' ^ ' 

where n is the density, g is 2, and 5 is a local moment. 

At high T increasing disorder, the number of on-site coupling of up and down spin 
electrons increases because electrons can occupy the low disordered sites and not occupy the 
strongly disordered sites. These coupled sites and unoccupied sites have no contribution on 
local moments, so x is niore suppressed with increasing disorder at high T. The plot of C 
justify this description. These are increased with increasing disorder as expected. On the 
other hand, at low T, we consider the Pauli paramagnetic susceptibility for free electrons: 

X''^^'' = filDie^), (15) 

where hb is the Bohr magneton and ep is the DOS at the Fermi level. For low disorder, there 
is a spin gap because DOS(eir) = 0. Increasing disorder, x is enhanced due to the increase 
of available states. Increasing disorder further (e.g. Vd = 4: or Vd = 5), this enhancement 
disappears because of the suppressed local moments at strong disordered sites. 

V. CONCLUSION 

In summary, we have calculated R,e[a{u!,T)], the frequency and temperature dependence 
of conductivity, for the two-dimensional Hubbard model with disorder at half filling using an 
inhomogeneous self-consistent Hartree-Fock numerical method. There is a gap in Re[(T(u;)] 
at low disorder. The characteristic temperature dependence of the conductivity in the Mott 
insulator shows -^ > for low T and -^ < at high T arising primarily from a carrier 
dominated regime at low T to a scattering dominated regime at high T. Upon increasing 
disorder, the spectral weight around the Mott gap is shifted to low frequencies and adc grows 
with disorder. (Xdc is enhanced at intermediate disorder with (Xdc/dT ^ for low temperature. 
However, adc is zero in the Mott insulator and in the strongly localized regime. 
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Moreover, we have calculated x(^ = 0,u;,T), uniform spin susceptibility, and Xs(^ = 
(vr, 7r),a;,T), staggered spin susceptibility. From the results of uniform static spin suscepti- 
bility, we considered the spin formation in each Region. 

We have investigated the effect of an inhomogeneous magnetic field, hi = — f^(cLcj|) 
considering the case with /ij = at all sites. We could have not seen any significant difference 
between two cases, hi = for all sites and hi ^ 0. Thus, the inhomogeneous field can only 
select the z direction on the 2D system and physical quantities that do not depend on z 
direction are not depending on hi. 

It may be interesting to contrast the behavior of the conductivity for the positive-U and 
negative-U disordered Hubbard model. Moreover, away from half filling is interested. 
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